Raman Scattering Study of the Lattice Dynamics of Superconducting LiFeAs 
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We report an investigation of the lattice dynamical properties of LiFeAs using inelastic light 
scattering. Five out of the six expected phonon modes are observed. The temperature evolution of 
their frequencies and linewidths is in good agreement with an anharmonic-decay model. We find 
no evidence for substantial electron-phonon coupling, and no superconductivity-induced phonon 
anomalies. 



PACS numbers: 74.70.Xa, 74.25.nd, 74.25.Kc 
I. INTRODUCTION 



Following the recent discovery of superconductivity in 
F-doped LaFeAsO with superconducting transition tem- 
perature T c ~ 26 K, 1 several families of Fe As-based su- 
perconductors including REFeAs(Oi- x F x ) (1 Ill-family, 
RE = rare earth), MFe 2 As 2 (122-family, M = Ba, Ca, 
Sr, K, Cs ...), MFeAs (Ill-family, M = Li, Na) 2 4 have 
been found and investigated. All these compounds share 
a similar tetragonal structure based on FeAs layers. In 
the 1111 and 122 compounds, superconductivity emerges 
when the structural and spin density wave transitions in 
the stoichiometric parent compounds are suppressed by 
chemical doping or pressure. LiFeAs deserves special at- 
tention, because it shows neither structural nor magnetic 
phase transitions, and superconductivity with T c ~ 18 K 
is present at ambient pressure without any dopingPHU 
The origin of superconductivity in this compound is still 
controversial. 



The presence of weak local moments^ and normal- 
state antiferromagnetic (AF) fluctuations jSEU combined 
with prediction of weak electron-phonon couplin^^ seem 
to support an unconventional pairing mechanism, as in 
other FeAs compounds. On the other hand, London 
penetration depth and small-angle neutron scattering™ 
angle-resolved photoemission spectroscopy (ARPES)p^ 
microwave surface impedance^ and several other mea- 
surements^ indicate a superconducting gap ratio 2A ~ 
4 ksT c closer to the weak-coupling limit than other 
Fe-based superconductors. In addition, ARPES mea- 
surements have been interpreted as evidence of strong 
electron-phonon coupl ing possibly enhanced by the AF 
spin fluctuations ! 19 * 20 ^ Taken together, these data rather 
suggest conventional electron-phonon driven pairing in 
this compound. The present study is motivated by the 
absence of experimental investigations of the lattice dy- 
namics that would allow a direct determination of the 
electron-phonon coupling strength in LiFeAs. 



II. EXPERIMENTAL DETAILS 

Single crystals of LiFeAs with T c ~ 18 K were grown in 
a sealed tungsten crucible using the Bridgman method.^ 
Since LiFeAs crystals are extremely air sensitive, they 
were cleaved and mounted on the cold finger of a helium- 
flow cryostat in a glove box under Ar atmosphere. The 
Raman spectra were taken in backscattering geometry 
in a JobinYvon LabRam 1800 single-grating spectrome- 
ter equipped with a razor-edge filter and a Peltier-cooled 
CCD camera. We used a linearly polarized He + /Ne + 
mixed-gas laser with a wavelength of 632.817 nm for ex- 
citation. The laser beam was focused through a 50 x 
microscope objective to a ~ 5 /im diameter spot on the 
sample surface. The power of the incident laser was kept 
less than 2 mW to avoid laser-induced heating. In or- 
der to determine the precise phonon frequencies at every 
temperature, Ne emission lines were recorded between all 
measurements. For data analysis, all phonon peaks were 
fitted by Lorentzian profiles, convoluted with the spec- 
trometer resolution function (a Gaussian line of 2 cm -1 
full width at half maximum (FWHM)). 



III. EXPERIMENTAL RESULTS 

A. Mode Assignment 

The lattice symmetry of LiFeAs is described by the 
space group P4/nmm (-Dj/J, with the Li, Fe and As 
atoms located at the 2c, 2b and 2c Wyckoff positions, 
respectively. From group symmetry analysis, 15 zone- 
center optical phonons are expected, among which 2Ai p , 
1 Bi g and 3E^ modes are Raman active! 2 ^ The corre- 
sponding atomic displacement patterns are sketched in 
Fig. [l]-a and the Raman selection rules as well as th e fre- 
quencies calculated within density functional theor;^^ 
are listed in Table 1. In Fig. [l]-b, we show the Raman 
spectra measured at room temperature for several scat- 
tering geometries with incident light wave vectors along 
the c-axis (z(xy)z, z(x'y')z, z(x f x f )z, z(xx)z configura- 
tions in Porto notation) or along the y axis (y(zz)y and 
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FIG. 1: (Color online) a) Atomic displacement patterns of the Raman-active optical modes of LiFeAs. b) Room temperature 
Raman spectra in z(xy)z, z(x'y')z, z(x'x')z, z(xx)z, y(xz)y, and y(zz)y configurations (note that the spectra in z(x'y')z and 
z(x'x')z configurations have been recorded in a lower-resolution mode to maximize the signal). Spectra have been shifted 
vertically for clarity. Starting from the left, the peaks are assigned to Ai^(As), Bi^(Fe), E p (l), E p (2), and Ai^(Li) vibrations, 
c) Aip (As) and Bi p (Fe) phonons for selected temperatures (room temperature, T ~ T c , and the base temperature of our 
cryostat). Black squares are the data, the red line is the result of a fit following the procedure described in the text. The 
spectra have been shifted vertically for clarity, d) Same plot for E^(l), E p (2), and Ai^(Li) phonons. 



y(xz)y). As expected, no phonon modes are observed in 
the z(xy)z geometry which selects the B2 g symmetry. 

The measurements in the z(x f y f )z, z(x f x f )z, z(xx)z 
configurations select phonons in the B^, Ai g and 
Ai^+Bi^ symmetries, respectively, which allows us to un- 
ambiguously assign the modes at 181 cm -1 and 223 cm -1 
to the Aip (As) and Bi^(Fe) c-axis polarized vibrations of 
the FeAs planes. These frequencies are in better agree- 
ment with the calculations of Ref. 13, which are based 
on the experimental lattice constants and atomic coordi- 
nates, than with those of Ref. 19, which were obtained 
using the relaxed ones. The mode frequencies are also 
very close to those of analogous phonons in the i2SplH26] 
or 1111 compounds pMS] We note that the second Ai g 
mode, involving mainly motions of the Li atoms along 
the c-axis, is not observed here. 

For measurements with light polarization in the ac 
plane, we observe three modes at 287 cm -1 , 297 cm -1 



and 328 cm -1 , as well as two other features at high en- 
ergy around 585 and 620 cm -1 . In the y(xz)y geometry, 
only the in-plane modes are expected. The 328 cm -1 
phonon is very weak in this configuration, but becomes 
much more intense than the two others upon switching 
to the y(zz)y configuration, and should therefore be at- 
tributed to the second A± g mode rather than to an in- 
plane E g phonon. The proximity of this mode frequency 
with the calculation for the Ai^(Li) mode further con- 
firms this assignment. The two remaining modes at 287 
cm -1 and 297 cm -1 are finally attributed to E^ phonons, 
despite significant disagreement with the calculated fre- 
quencies (240 (225) and 290 (233) cm" 1 from Ref. [13] 
and \19\ respectively). We note that the selection rules 
are not perfectly respected here as, for instance, the E^ 
modes are still visible in the y(zz)y geometry. The ac 
plane is not a good cleavage plane, giving rise to a rough 
surface hard to align accurately in the glove box. As 
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TABLE I: Calculated Raman-active phonon frequencies and selection rules from Refs. [13] and [19] and comparison to our 
experimental results (see text for the definition of the parameters). 



Mode 


Polarization 


Selection Rule 


Calc. Freq. (cm 1 ) 


Calc. Freq. (cm 1 ) 


Exp. fitting Parameters (this 


work) 


(atom) 






from Ref . [13] 


from Ref. [T9l 


ojq (cm 1 ) 


C (cm" 1 ) 


To (cm" 1 ) 


T (cm -1 ) 


E g (As) 


in-plane 


xz 


121 cm" 1 


173.7 










Ai g (As) 


c-axis 


xx, zz, x'x' 


188 


183.3 


185.1 


1 


0.24 


1.3 


Bi 9 (Fe) 


c-axis 


x y , xx 


225 


207.5 


237.8 


4 


0.07 


1.6 


E g (Fe) 


in-plane 


xz 


240 


224.7 


289.2 


0.9 





1.1 


E g (Li) 


in-plane 


xz 


294 


233 


308.2 


3.7 
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338.2 
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FIG. 2: (Color online) Comparison of the temperature depen- 
dence of the frequencies of the high-energy modes at 585 and 
620 cm -1 with the sums of the frequencies of modes observed 
in y(zz)y and y(xz)y polarizations. 



a consequence we always get a slight misalignment in 
the crystal that causes the observed mode leakages. The 
lowest-energy mode predicted by these calculations could 
not be observed in this study. It is worth noting that 
the two Aip modes are not visible in the same scatter- 
ing geometry. In z(xx)z or z(x"x')z configurations, the 
phonon intensity is proportional (2a xx ) 2 , and to (a zz ) 2 
in the y(zz)y configuration, where a zz and a xx are the 
diagonal elements of the A lg Raman tensor. Our results 
imply that the Raman tensor components for the two Ai g 
phonons are different: we have indeed a zz ~ << a xx 
for the Aig As mode while we have a xx ~ << a zz for 
the Aig Li mode. Such anisotropy of the Raman tensor 
has for instance been reported in the case of SrFe2As2p' 
where a xx ~ << a zz has been observed for the Ai g 
mode. 



The energy of the two features observed at 585 and 620 
cm -1 is much higher than the highest calculated p hono n 
energy in this system (which is about 300 cm -1 UJiLS]). 
They can therefore not be attributed to single-phonon 
modes, but rather appear to be due to two-phonon scat- 
tering. This is further confirmed by the observation that 
585 cm -1 is almost exactly equal to the sum of the fre- 
quencies of the two Eg modes at all temperatures, while 
the frequency of the second feature is always located be- 
tween the sums of the energy of the Ai p (Li) phonon with 
each of the two modes (Fig. |2j. We note that the 
large width of this feature compared to the one at 585 
cm -1 (~ 35 cm -1 vs. ~ 3.5 cm -1 ) indeed suggests that 
it can be further decomposed into several features that 
are not individually resolved. This is not clear why these 
features are observed only when the light is polarized in 
the ac plane, and further investigation are required to 
understand this point. 

As seen in Figs. [I}c and d, all of these phonons are 
found to harden and sharpen as the temperature de- 
creases. Note that the intensity drop seen between 20 
and 5 K in Fig. [l]-d was not reproducible and shall 
not be considered as a superconductivity-induced effect. 
Phonon intensities will not be further considered. Fig- 
ure [3] shows the temperature dependence of their posi- 
tions and linewidths. The absence of phonon anomalies 
at any temperature (such as the splitting of the E^ modes 
predicted^ and observed^ at the structural phase tran- 
sition in BaFe2As2) confirms that LiFeAs does not un- 
dergo any structural or magnetic phase transitions. Fig- 
ure [3] also shows the results of fits of the experimental 
data points to a simple expression for the temperature 
dependence of the frequency and linewidth of an optical 
phonon resulting from symmetric anharmonic decay, i.e. 
decay into two acoustic modes with identical frequencies 
and opposite momenta! 32 * 33 * 



(T) = oj -C 



T ph (T) = T + T 



1 



1 



(1) 



(2) 



e 2k B T _ l_ 

where C and T are positive constants, ujq is the bare 
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FIG. 3: (Color online) Temperature dependence of frequen- 
cies (panels a-e) and linewidths (panels f-j) of the five observed 
Raman-active modes. The pink dashed line marks T c , and the 
green line is the result of a conventional phonon anharmonic 
model (see text). 

phonon frequency, and To a residual (temperature in- 
dependent) linewidth originating from sample imperfec- 
tions or electron-phonon interactions. We find excellent 
agreement with this expression. The fitting parameters 
are summarized in Table 2. 

In all cases, the residual width r is found to be van- 
ishingly small, which confirms the high quality of our 
crystals. The temperature dependent coefficient T is al- 
ways much larger than the residual width, indicating that 
lattice anharmonicity is the principal source of decay for 



the phonons in LiFeAs, in clear contrast with the recent 
observations of Litvinchuk et al. in Pr x Cai_ ;E Fe2As2!^' 
Three kinks were recently observed around 15, 30 and 
44 meV (121, 242 and 355 cm -1 ) in energy distributions 
curves measured by ARPES on LiFeAs, and attributed 
to electron-phonon coupling.^ Our present data cannot 
completely rule out the possibility of such a coupling, 
but excludes it to occur at the zone center where our 
measurements were carried out. To confirm the electron- 
phonon coupling origin of this kink, a study of the phon on 
dispersion is this compound is therefore required- 34 * 35 *. 

No changes in any of the phonon frequencies or 
linewidths are observed at the superconducting transi- 
tion, as previously reported for K-doped BaFe2As2r^and 
F-doped NdFeAsO.^ This is consistent with the rela- 
tively small amplitude of the superconducting gap (2A ~ 
4.0 ksTc ~ 50 cm -1 compared to the phonon frequen- 
cies reported hereP To our knowledge, the only Fe As- 
based compound in which superconductivity-induced 
changes in the Raman phonons are clearly observed is 
Pr^Cai-^FesAssP 



IV. CONCLUSIONS 

In conclusion, we have reported a lattice dynamical 
study of superconducting LiFeAs. We observed five 
modes out the six expected. Their frequencies are in good 
overall agreement with values predicted by density func- 
tional calculations, and their temperature dependence is 
well described by an anharmonic-decay model. We ob- 
served neither clear anomalies associated with the super- 
conducting transition nor any evidence for substantial 
electron-phonon coupling. 
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